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Shor's factoring algorithm illustrates the potential power of quantum computation. 
Here, we present and numerically investigate a proposal for a compiled version of 
such an algorithm based on a quantum- wire network exploiting the potential of fully 
coherent electron transport assisted by the surface acoustic waves. Specifically, a 
non standard approach is used to implement, in a simple form, the quantum circuits 
of the modular exponentiation execution for the simplest instance of the Shor's 
algorithm, that is the factorization of iV=15. The numerical procedure is based 
on a time-dependent solution of the multi-particle Schrodinger equation. The near- 
ideal algorithm performance and the large estimated fidelity indicate the efficiency 
of the protocol implemented, which also results to be almost insensitive to small 
destabilizing effects during quantum computation. 
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I. INTRODUCTION 

Quantum computers can efficiently solve some problems that are unaffordable on clas- 
sical computers, and processing the information encoded in quantum systems results to be 
extremely powerful for particular tasks. Specifically, quantum mechanics effects like entan- 



*Elcctronic address: fabrizio.buscemi@unimore.it 



2 



glement and wavefunction superposition turn out to be fundamental building blocks, and 
allow for the quantum computational speedup over classical computation. The Shor's algo- 



rithm 



LH4j has undoubtedly been widely investigated among those illustrating the power of 
quantum computation. As a matter of fact, it plays a key role in cryptographic protocols, 
because it allows one to factorize a composite number with a computational time that is a 
polynomial function instead of an exponential function of the number itself. 

The practical implementation of the Shor's algorithm represents a challenge for quan- 
tum information science. Two possible physical architectures have been proposed: nuclear 
magnetic resonance (NMR) 0,0] and photonic 7-^] systems, even though some open ques- 
tions exist in both cases. While in NMR it is difficult to prepare the qubits in pure states 
and control their coherent evolution, thus leading to a controversial quantum nature of the 
experiments, photonic systems cannot be scaled to a larger number of qubits due to their 
size and stability limitations. Nevertheless, a recent experimental demonstration of Shor's 
algorithm was obtained by means of optical waveguides integrated on silica-on-silicon chips. 
Even if the efficiency of the single photon source and detectors still does not appear to be 
very good {9], the suggested architecture is promising for the implementation of large-scale 
quantum circuits on many qubits. 

No evidence of compiled version of quantum factoring algorithm using electron qubits 
has been achieved so far. The approach of using charge carriers in solid-state systems is 
very appealing because it does not allow only to overcome the scalability problem, but it 
also provides a valid guideline for the design of devices easily integrable in the traditional 
electronic circuitry. Specifically, the possibility to implement Shor's quantum factoring 
algorithm on an electronic chip would certainly represent an essential test to verify the 
potential of quantum cryptography on everyday life. 

In this paper, we propose and numerically simulate a compiled version of the Shor's algo- 
rithm. Electronic quantum logic gates in one-dimensional (ID) semiconductor channels are 
used to realize the necessary processes and to produce multi-particle entanglement and mul- 
tipath interference. In particular, we have considered the fully coherent Surface-Acoustic- 
Wave (SAW) assisted electron transport in couples of GaAs quantum wires lOj, with the 
qubit defined by the localization of a single carrier in one of the coupled channels 11]. 

Quantum-wire systems have been shown to be suitable to produce bipartite entangled 
states [12] and to perform quantum teleportation Here, the numerical implementation 



3 



of the quantum factoring algorithm results to be much demanding in comparison with the 



previous works 



12 



13( , due to the higher number of the simulated quantum logic operations 



over many qubits. Specifically, we design the quantum circuits of the modular exponentiation 
execution for the easiest meaningful instance of Shor's algorithm, that is the factorization 
of A=15 for two different co-primes (7=11 and C=2 (defined in Sec. HI!), corresponding to 
the period r=2 and r=4, respectively. The circuit performing the modular exponentiation 
function is brought to a form different from the one given in literature l|, 4]. This procedure 
allows one to move on toward simpler networks of electron quantum gates, and aims at 
future research leading to a scalable full-realization of the Shor's algorithm in quantum- 
wire devices. In our implementation the inverse quantum Fourier transformation (QFT) is 
not present since it is not necessary for any order-2 z circuit (with I 6 N), as shown in the 
literature [7] . For sake of completeness, a description of the circuit realizing the inverse QFT 
in a quantum-wire network has been given elsewhere jl4]. 

This paper is organized as follows. In Sec. [TT] we illustrate the theoretical features of 
Shor's algorithm, while the description of the physical implementation in a quantum-wire 
device and the discussion of the numerical approach adopted are given in Sec. IIHI In Sec. IIVI 
we show the results obtained from the numerical simulations of two quantum circuits for 
the factorization of iV=15 corresponding two different parameters choice. Comments on the 
results and final remarks are drawn in Sec. |VJ 



II. SHOR'S ALGORITHM 



The strategy to find a nontrivial prime factor of the positive integer A is described in 
the following. A random co-prime C is chosen, i.e. A and C have not common factors. 
Euler's theorem states that exists an integer r such that C r = 1 mod A (that is, C r — 1 is an 
integer number multiple of A), with 1 < r < A. The number r is called order of C mod A. 
Provided that the latter is even, then it follows that C r -1 = (C r/2 -l){C r/2 + l) = mod A 
and this implies that A is a divisor of the product (C r / 2 — l)(C r / 2 + 1). Assuming that 
(C r / 2 7^—1 mod A), it follows that A must have a common factor with both (C r / 2 ± 1). 
Therefore, this implies that the factors of A are given by the greatest common divisor of A 
and (C r / 2 ±1), which can be efficiently computed by means of Euclid's classical algorithm. It 
is worth noting that in order to guarantee the algorithm validity the two conditions stating 
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that r is even and {C T I 2 7^ — 1 mod N) must be satisfied. These conditions are met with 
high probability for iV odd, except in the case where N is a prime power (N = p a with p 
prime). Thus, the smallest composite integer N that can be successfully factored by Shor's 
method is iV = 15. If N is an even number or a prime power, other classical methods should 
fruitfully be applied for the factorization instead of the Shor's method. 

The Shor's algorithm needs quantum computation only in the evaluation of the order 
r. The quantum order-finding routine commonly uses two registers of qubits the 
argument register with n = 2 hi2 N qubits, and the function-register with m = hi2 N qubits. 
Its implementation can be separated into three distinct steps. The first one is the register 
initialization corresponding to 

| )®«|0)® m 1—). -L (|0) + \l)f n \0)® m - l \l) = 
v 2 n 



2 n — 1 

-L J2 Wlo) " 1 - 1 !!}, (i) 

where the argument register is prepared by means of the Hadamard transformations in an 
equal superposition of all n-qubit computational basis |0(l) a;i )|0(l) a ; 2 ) • • • 10(1)^) • • • |0(l) In ). 
In the second step, also known as modular exponentiation, the function C x mod N is imple- 
mented on the function register, while the argument register remains in x. The global state 
is thus given by 

1 2 n -l 

-=5>>|CWIiV>. (2) 

The state of Eq. t$Z§ is highly entangled and exhibits the so-called "massive parallelism" , i.e. 
the execution entangles in parallel all the 2™ input values with the corresponding values of 
C x 'mod N, although the algorithm has run only once 0, [ijj] . Finally, the inverse QFT is 
applied on the argument register yielding the state 

2 n -12 n -l 

E eW2 >>i^ modiV >> ( 3 ) 

where, due to the interference, only the terms \z) with 

z = oTjr (4) 



have a significative amplitude. In the above, a is a random integer ranging from to r-1. 
Thus, if one performs measurements on the outcome of the argument register, he would 
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get a2 n /r for some a, and the order r can be deduced after the classical procedure with 

n 

probability greater than 1/2 (see Ref. [3j). 

The modular exponentiation, that is the evaluation of C^mod N for 2 n values of x in 
parallel, is the most demanding part of the algorithm. This can be performed by using the 
identity: 

2 (n-l) 

+ ■ • • + x\2 l + xo2°, where Xk are the binary digits of x. From this, 

it follows that 

C x mod N = C 2 ^'^ 1 ■ ■ ■ C 2xi C x °mod N = 
C 2 '"" 1 ' 1 - 1 • • ■[C 2xi [C Xo mod iVjmod N) ■ ■ ■ mod N]. (5) 

This means that we first multiply 1 by Cmod N, if and only if xq = 1; then we multi- 
ply the result by C 2 mod N if and only if x\ — 1 and so forth, until we finally multiply 
by C 2< ~" 15 mod N if and only if = I. Therefore, the modular exponentiation consists 
of n serial multiplications modulo N, each of them controlled by the qubit x^. The fac- 
tors C, C 2 , ■ ■ ■ , C 2( " ^mod iV can be found efficiently on a classical computer by repeated 
squaring. 



A. N=15 



As noted in the previous section, the Shor's factorization algorithm fails if N is even or 
a prime power, and the smallest composite integer that can be successfully factored by 
means of Shor's method is A^=15. Even if N is small, this compiled version of the Shor's 
algorithm displays a great potential for a future realization of large-scale quantum algorithm. 

Being A^=15, the minimum size of the function and argument registers must be 
m=ln 2 [15]=4 and n=2m=8, respectively. The algorithm would then require at least 12 
qubits. However, the following comments allow us to reduce the number of qubits neces- 
sary for the purpose of a proof-of-principle demonstration. A co-prime C with 15 is one 
element of the set 2, 4, 7, 8, 11, 13, and 14. As shown in Table [H it comes from repeated 
squaring that C 4 mod 15 = 1 for all valid C. In turn, this implies that only two bits Xq 
and X\ are needed for the controlled multiplications. As a consequence, the multiplications 
by C 4 ,C 8 ,--- are trivial, and all the multiplications, except the ones by C and C 2 , can 
be left out. For C = 4, 11, 14, C 2 mod 15 = 1 and only the first bit xq is relevant. These 
considerations account for a reduction of the size of the argument register, which can finally 
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TABLE I: The table displays C^mod 15 for all C < 15 co-prime with 15 and for values of x which 
are power of two. Note that C^mod 15 = 1 for all valid C. 

be constituted by no more than two qubits (n—2). Adding this latter result to the four 
qubits of the function register, only six qubits are needed instead of twelve, as previously 
found. 

Shor's factorization algorithm for the number 15 turns out to be particularly simple 
because it does not require the implementation of the inverse QFT in the quantum circuit. 
As shown in the literature jj], the latter is not necessary for any circuit of order 2 l and It can 
be replaced by a classical processing which also inverts the order of the computed quantum 
bits of the argument register. 

In this work we implement the Shor's quantum factoring algorithm and check it against 
two co-primes: C=ll and C=2, that are representative parameters for the system at hand. 



1. C=ll 



This parameter choice represent an "easy case" since the modular exponentiation can be 
simplified to the multiplication of the initial function register state, |y) = |0 y3 0j /2 0j /1 lj /0 )=l, 
by C=ll controlled only by Xq js]. 

In the left panel of the Fig. [I] a compiled version of the quantum circuit for (7=11, using 
the inverse QFT, is displayed. At first, both registers are initialized: each qubit of the 
argument is prepared by Hadamard gates in a superposition of and 1, and the function 
register state is set to \y)—l, so that the global state |$c=n) °f the system is 



|$C7=ll) = 7j (l°*Ao> + lOxxlxo) + |lxA ) + IWso)) 

X l y 3 0y 2 ly ) . 



(6) 
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Then, the modular exponentiation is performed: the controlled multiplication of 1 by 11 
is equivalent to the controlled addition of 10 to 1. The latter is implemented in the quantum 
circuit by two controlled-NOT (CNOT) gates: one between xq and y\ and one between xq 
and y%. It is worth noting that the qubits yo and y^ evolve trivially during computation. 
Thus the state of the system takes the form 

3 

|$ c=11 ) = ^ \x)\ll x mod 15) = 

x=0 

2 X M^a;i0a!o)l^»8^ya^»i-'-l/o) \^xi^-x )\^-y s O y2 l yi ly ) 

+ |lxiO XO ) |0y 3 0y 2 yi ly Q ) + |l Il l :ro ) |ly s 0j /2 l ?/l lj /o )^ . (7) 

This means that a Greenberger-Horne-Zeilinger (GHZ) entangled state 

l*> = -^(|0^0«iO»i> + lW«,l«i» (8) 

is created between qubit xq of the argument register and qubits y% and y\ of the function 
registers. 

The final step is represented by the inverse QFT. The right panel of the Fig. [T] 
shows a compiled version of the quantum circuit without using the inverse QFT. Please 
note that in this case also qubit x\ is redundant: the corresponding Hadamard gate 
results to be unnecessary and does not need to be implemented. Here, the ini- 
tial state is 4g (lO^O^) + lO^l^)) lO^O^O^l^) and the modular exponentiation yields 

"^TJ (|0a;i0a; o )|0j /3 0y 2 0j /1 ly o ) + lO^l^g) | ly 3 0y 2 ^-yi^-yo))- 

The outcomes of the measurement on the inverted argument qubits Xq and x\ give then 
00 or 10 with equal probability. Once this result is known, one can obtain the order r of C 
mod iV from Eq. fll]). While the output 00 corresponds to a failure, the output 10 allows 
one to determine the period r = 2 2 /2 = 2 and represents a successful implementation of the 
Shor's algorithm. 

2. C=2 

Since the number of gates needed to perform the modular exponentiation is greater than 



the case of C=ll p|, the choice of C=2 represents a "difficult case". In fact, the modular 
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FIG. 1: (Color online). Left panel: Outline of the quantum circuit for quantum factorization of 15 



for (7=11, using the inverse QFT. 



B] Right panel: Outline of the quantum circuit for quantum 



factorization of 15 for C=ll, not using the inverse QFT. 

exponentiation is given by the multiplication of \y)=l by 2 controlled by xq and by the multi- 
plication of the obtained result by 4 controlled by x%. The left panel of Fig. [2] shows the quan- 
tum circuit for the case at hand. The network for the modular exponentiation is composed 
of two CNOT followed by two controlled SWAP (CSWAP) gates; the first two correspond to 
the addition of 1 to |y)=l controlled by xq, while the CSWAP gates multiply the result by 4 
controlled by x\. The modular exponentiation leads to the state: J2 x =o |^)|2 a: mod 15) = | x 

( | 0^! O^) \0y 3 0y 2 yi ly ) + | 0^ 1 Xq ) \ 0y 3 0y 2 1 y y 0y Q ) + | 1 Xj Xq ) | y 3 0y 2 1 0y Q ) + | 1 Xl ^X ) | 1 y 3 ^mPtll ^VO ) ) ' 

Nevertheless, a different compilation of the q nantnm circnit can be realized Q, and it 
is reported in the right panel of Fig. [2j By means of the latter, it is possible to evaluate 
ln 2 [2 x 'mod 15] in the function register in place of 2 x mod 15, thus reducing the number 
of function qubits from ln 2 [15]=4 to ln 2 {ln 2 [15]}=2. This compilation maintains all the 
features of the algorithm originally proposed 3] , and still does not make use of the inverse 
QFT, as in the previous case. Following this scheme, the initialization of the system leads 
to the state 



|$C=2> 



X1 XO ) + \0 X1 1 XO ) 



+ |lxiO Xo ) + ll-cilaso) ) I Oyi Oyo ) 



(9) 



meaning that the argument register is kept in the usual equally-weighted coherent superpo- 
sition of all possible arguments, while the initial function register state is \y)=0. If we apply 
the procedure described in Eq. ((Sj) to evaluate ln 2 [2 z mod 15], it can be easily shown that 
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the modular exponentiation reduces to the sum of ln 2 [2 1 mod 15] = 1 to \y)=0 controlled by 
xq, and of ln 2 [2 2 mod 15]=2 to the obtained result controlled by X\. These operations are 
implemented in the quantum circuit reported in the right panel of Fig. [2] by means of two 
CNOT gates: one between xq and yo and another between x\ and y±. It is worth noting 
that in this case the algorithm is very simple since it consists of only two networks of gates 
acting on independent qubit pairs. After modular exponentiation the state of the system 
takes the form 

|$C= 2 ) - U\0x 1 0x 0y 1 0y ) + 1 Xl l XQ yi l yo ) 
+ | l^i 0x lj/x O^q ) + ll^l^gl-j^lyo) 

= \(\0x 1 in ) + IWwH (\0x %o) + IWwA ( 10 ) 

that is the product of two entangled Bell pairs, thus confirming the manifestation of entan- 
glement between the two registers of the algorithm. The inverse QFT is not necessary and 
it can be replaced by its classical counterpart, which also swaps the output quantum bit of 
the argument register. The two-bit outputs for the case under investigations are: 00, 01,10, 
and 11. The second and the fourth outcomes allow the evaluation of the order r=4, which 
efficiently yields the factors 3 and 5 via Euclid's classical algorithm; the first one corresponds 
to a failure mode and, lastly, the third one leads to trivial factors. 

III. THE PHYSICAL IMPLEMENTATION AND THE NUMERICAL 

APPROACH 

Here, we describe the implementation of the Shor's algorithm in a specific semiconductor 
nanostructure. It consists of a number of couples of GaAs quantum wires where surface 
acoustic waves (SAWs), i.e. sinusoidal piezoelectric potential, propagate and trap charged 



carriers into their moving minima, letting one particle fill in each minimum |16|. The so- 
called flying qubits are realized by means of the states |0) and |1), encoded through the 
localization of a single electron in one of the two ID channels [U|. Here, the SAWs are used 
to inject and drive the electron thanks to their efficiency in preventing the natural spatial 
spread of the wavefunction 1t| and in making the carriers more immune to the decohering 
effects 18]. Moreover, in this investigation the carrier transport is assumed to be fully 
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FIG. 2: (Color online). Left panel: Outline of the quantum circuit for quantum factorization of 15 
for C=2, using the inverse QFT and evaluating 2 2 mod 15 in the function register. Right panel: 
Outline of the quantum circuit for quantum factorization of 15 for C=2, not using the inverse 
QFT and evaluating ln2[2 a hiod 15] in the function register. Both the circuits correspond to the 



ones examined in Ref. 



coherent. 

As shown in the literature 



19 



20J, such a system is able to provide the universal set of 
gates useful to realize any quantum computational network. Specifically the basic building 
blocks are R x (9), i?o(i)(0), an d T(j) \2l\ . The former two gates implement one-qubit logical 
operations, whereas the latter is a two-qubit gate. 

R x (8) acts as an electronic beam splitter and can be materialized through a coupling 
window between the two wires of the qubit [221 ]. Its matrix representation on the basis 

{I ))! 1 )} is § iven b y 



RM 



cos 2 i sin | 



i sin 2 cos 2 



(11) 



Ron) ((f)) is an electronic phase shifter obtained by inserting a potential barrier in the 
wire 0(1), thus inducing a delay phase in the propagation of wavefunction. Its action is 
described in one-qubit basis by 



Ro(<f>) 




and 



RM) 



1 o 




(12) 



T(7) is a conditional phase gate exploiting the Coulomb interaction between two electrons. 
It consists of a region in which the carriers propagate along two different wires close enough 
to give rise to an effective interaction able to delay both particles. The matrix representation 
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of T( 7 ) in the two-qubit basis {|00),|01), |10), |11)} is: 

( I o\ 
10 

T{l) 



(13) 




\0 1 / 

The phases 9, 0, and 7 of the above quantum gates depend upon the physical and geometrical 
parameters of the systems such as velocity, amplitude and wavelength of the SAW potential, 
strength of the electron-electron interaction, coupling window length, and shape of the 
potential barrier. In order to perform any transformation of the many-qubit state, an 
appropriate tuning of the above parameters in a given network of R x (9), i? O (i)(0); an d 
T( 7 ) gates is required. 

In the compiled versions of the Shor's algorithm illustrated in the right panels of Figures [D 
and EJ corresponding to C=ll and C=2, respectively, the two logical operations involved 
are only the Hadarmad H and CNOT gates. In terms of R x (9), Ro(i)(<P), and T(j), these 
operations can be reworked as [14] : 




(14) 



and 



CNOT = R, 



(2) 



^4,(|) r „, w « ? . ( |) 4 , ( |) fl a» w ._L 



( I o\ 

10 
1 
\00 1 0/ 



(15) 



The numerical implementation of the quantum circuits described in the previous section 
results extremely challenging due to the large number of basic building blocks needed to 
realize the H and CNOT gates. Furthermore, the experimental realization of such devices 
would also certainly meet serious obstacles stemming from the difficulty of preserving the 
coherent evolution of the qubits during a number of logical operations, from decoherence 
effects due to interactions with the environment, as well as from possible structural defects 
induced by the processes of fabrication and tuning of each quantum gate. The theoretical 
and experimental implementation of two-qubit quantum circuits by means of a minimum 
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amount of operations has been widely discussed in literature 23. |24J| . Different protocols 
have been proposed ranging from the special perfect entanglers 23[ to the use of a given 
tunable entangling interaction [24j. Here, we propose in the followings a scheme suitable 
to perform the quantum factoring algorithm in devices formed by semiconductor quantum 
wires with a minimal number of the fundamental gates R x (6), -Ro(i)(0) ; an d T ('-/). The 
proposed implementation satisfies the main requirements of the Shor's algorithm as originally 
formulated |lt|3|. In fact, "massive parallelism" is maintained, since entanglement is created 
between argument and function registers, and the binary output of the argument qubits are 
unchanged, as it will be shown below. In the left and right panel of the Fig. |3]we report the 
quantum- wire networks implementing the circuits displayed in the right panels of Figures [U 
and El respectively. 

In the first case, corresponding to (7=11, the network implemented reads 

R ( : o) (£) (I) I*»*>(7r)J#0 (|) R(vs) (£) tfW ^tM^rM {v)B to) Q ; 

(16) 

where the superscripts of the quantum gates indicate which qubit they act on. The three- 
qubit output state is: 

^~ l^oOjraOyi) + i\^-x ^y3^yi)) • (17) 
For (7=2 the global logical transformation can be written as: 

R { : o) (£) (f ) (|) Rt ] (I) T^^T^^R^ (I) Hjn) (I) (18) 

and, after being applied to the input state \0 Xl Xo yi 0y o ) , it yields 

2 — l^i^J/i)^ ^|0x o 0j/o) — I la;o lyo ) ^ • (19) 

In both cases, the degree of the entanglement created between the qubits of the argument 
and the function register is equal to the one of the standard formulation of Shor's quantum 



factoring algorithm 



|3|. While for (7=11 the output state is GHZ-like, i.e. the maximum 
amount of quantum correlations is built up among three qubits, when (7=2 the four qub it- 
state is the product of two entangled Bell pairs. Furthermore, the reduced density matrices 
of the argument qubits correspond to the ones calculated from the quantum states of Eqs. (JSJ) 
and ( flOl) . In turn, this implies that the binary output of the algorithm, i.e. the outcome 
measurements on the argument register, be identical. 
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FIG. 3: (Color online). Left panel: Sketch of the physical system used to factorize JV=15 with 
C=ll, corresponding to the quantum circuit displayed in the right panel of Fig. [TJ The beam- 
splitter R x of the qubit xq mimics the initialization procedure, that is the equal splitting of electron 
wavefunction of the qubit xq in the two wires. The next two sets of gates, first R X TR X acting 
on {xq^ui} and then R x RqTR\R x on {xo,y%}, play the role of two CNOT gates creating the 
maximum entanglement of the qubits {xo, yi, 2/3} in the modular exponentiation step. For sake 
of clarity the qubits yo and y2 evolving trivially during the computation have not been reported. 
Right panel: Sketch of the physical system used to factorize iV=15 with C=2, corresponding to 
the quantum circuit displayed in the right panel of Fig. [5J In the initialization procedure two gates 
R x act on the the couple of the argument register qubits {xq,xi} and split each of them in an 
equal superposition of the |0) and |1) states. The modular exponentiation consists of two network 
R X TR X , each operating onto {xo>2/o} an d {xi,yi} mimicking the action of two CNOT gates. Note 
that here, for brevity, the phases of the quantum gates involved in the networks and explicitly 
indicated in the Eqs. (fT6j) and (fTBj) . are omitted. 



The networks of the gates described in Eqs. (|T6|) and (|T8|) have been simulated by solving 
numerically the time-dependent Schrodinger equations for three and four electrons injected in 
the GaAs quantum wire devices of Fig. [31 While for C=2 the four-particle dynamics reduces 
to the time evolution of a pair of separable two-particle systems, when C—ll, a solution of 
the Schrodinger equation for the whole three-carrier wavefunction is required 25j. Anyway, 
in order to minimize the computational burden of both of the reported cases, a semi-ID 
approach has been used to investigate the time evolution of the system in place of a two- 
dimensional (2D) computational scheme. This approach was already introduced to simulate 
a teleportation protocol in a quantum wire device [ljj. According to this simplified scheme, 
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the direction y of the particles along the wires is fully included in the simulation (y being 
discretized with a point-grid resolution of Ay=lnm), while the two variables describing the 
longitudinal direction x and identifying the wire where the carriers are localized can only 
assume one of the two values or 1. 

Though it does not allow to simulate the gate R x (9), which would require a full 2D 
analysis, the numerical procedure here adopted makes it possible to move from a time de- 
pendent Schrodinger equation for a mult i- variable wavefunction $(X, Y, t) (seven unknowns 
for the device with C=ll and nine when C=2) to many coupled Schrodinger equations of 
the following kind: 

d , ,„ / h 2 



i^*x(Y, t) = [~^y + Vx(Y, t) $ X (Y, t), (20) 

where X = (x Xo ,x yo ,x Xl ,x yi , . . .) and Y = (y ao , y yo ,y Xl ,y yi , . . .). Specifically, when C=ll 
a system of eight coupled equations is obtained, while for C=2 two independent systems 
of four equations are found since the qubit pairs {xo,yo} and {x\,y\} are independent. In 
both the Crank-Nicholson finite difference scheme has been applied to get a nu- 

merical solution. The potential term Vx(Y, t) appearing in Eq. (120]) sums up the SAW time 
dependent potential, the Coulomb interaction between electrons, and the static potential 
profile. The simulations here presented make use of a sinusoidal potential mimicking a SAW 
of amplitude and wavelength equal to 20 meV and 200 nm, respectively, and propagating 
with the sound velocity v s = 3.3 x 10 3 m s _1 . Screening effects have been included in the 
Coulomb potential between the carriers by inserting an exponential damping term 27| with 
a Debye wave vector of 0.2 nm -1 . 

In order to numerically implement the networks described in Eqs. ([16]) and ([18]) . one must 
firstly to find the suitable geometrical parameters of the device for the gates i? O (i)(0) and 
T(^y) giving the required value tt of the phases <j) and 7. To this aim, we have performed 
a number of simulations testing different geometries for the phase shift and the conditional 
phase gate. As to the phase <fi is concerned, it depends on the height and the length of the 
potential barrier. The values of these parameters obtained from the optimization procedure 
are 2.82 meV and 8 nm, respectively, and correspond to a delay phase <fi of 0.927T, that is 
good enough for our purposes, as it will be shown in the next section. It is worth noting 
that the barrier height is significantly smaller than the amplitude of the SAW potential, this 
making the spatial spreading of the electron wavepacket negligible and letting it be entirely 
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transmitted through the barrier. The main geometrical parameters affecting the phase 7 of 
the conditional phase gate are the length of the coupling region and the distance between the 
coupled wires. Their optimal values used in the numerical implementation of the algorithm 
are 150 nm and 5 nm, respectively. They lead to a 7 value about 0.887T, which allows us 
to simulate satisfactorily the two-qubit logical operation of a CNOT-like gate. By applying 
the described geometry for T(^y), both the tunnelling effects between the two wires and the 
reflection phenomena in the coupling region are negligible. 

From a computational point of view, the numerical simulation of the T(^y) gate is more 
challenging than that of the phase shift i2o(i) (</>)■ While the latter involves a one-particle 
potential, the former exploits a two particle-interaction that builds up an amount of quantum 
correlations between the wire degrees of freedom of the particles, as expected. However, it 
creates also an undesired entanglement between the variables defining the position of the 
carriers along the wires. As a consequence, the evaluation of the effects of the controlled 
phase gate on the multi-particle wavefunction is a demanding task because it implies that a 
number of two particle simulations must be combined together to obtain the time-evolved 
state of the overall system. 

The R x (6) gate has not directly been simulated; nevertheless, its action has been taken 
into account by means of the transformation matrix of Eq. (TTT!) . validated by the results of 



appropriate 2D simulations 



20). 



IV. RESULTS 

According to the sectioning of the previous paragraphs, we firstly present the results 
obtained for C=ll, then those for C=2. 



A. C=ll 



Figure H] shows the density matrix describing the qubits of the argument and function 
register, {xo, yi, y%} at three different stages: input, initialization procedure and modular 
exponentiation. For sake of simplicity, we do not consider the x\, yo, and 2/2 qubits that 
evolve trivially during the computation. The knowledge of the joint state of both registers 
after modular exponentiation is essential for the estimation of the device performance. In 
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b) 




c) 



FIG. 4: (Color online). The density matrix of the qubits {£o> 2/1)2/3} evaluated at three different 
time steps: a) input, b) initialization procedure, and c) after the modular exponentiation. Note 
that these matrices have been obtained from the full density matrices of the electrons of the 
qubits {xo, 2/1,2/3} by integrating over the variables defining the position of three carriers along the 
channels. Here the moduli of the density matrix elements are plotted. 



particular, we find that the output quantum state corresponds to the GHZ-like entangled 
state of Eq. (TTT|) with a good approximation. For a more quantitative evaluation of the 
reliability of the algorithm, we have also calculated the fidelity F = (fy\p out \^), where p ou t 
is the output density matrix of the full system. The high value found, F=0.97, evidences 
the very good efficiency of the implementation. Such a result is certainly related to the fact 
that our numerical simulations have been performed by setting the device temperature at 
K, that is neglecting any effect of decoherence induced by the environment on the carrier 
transport properties. In particular, this means that the electron-phonon interaction have 
not been included in the simulations. These assumptions are physically sound when we take 
into account that the experimental investigations of the low-dimensional structures used as 



the basic blocks for our device are usually performed at very low temperatures [lfj, 122 ] . 
Moreover, one of the pros of our results certainly relies on the high fidelity value, which has 
not been obtained under ideal geometries for the i?o(l)(0) and T(j) gates. Once more, this 
off-ideality situation is close to the experimental conditions. 

The density matrix of the argument register after modular exponentiation is displayed 
in Fig. [5j Specifically, we show the reduced density matrix p{ XQ)X >}(y, y) of the electron of 
qubit xo, without the redundant qubit x\. The output of the quantum circuit is the logical 
state probability, that is the probability of finding the electron of the qubit xq in the wire 
or 1. This is described by the integral over y of the diagonal elements of p X0;X '(y, y). 
The off-diagonal terms are very small, thus proving that the argument register becomes a 
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full quantum statistical mixture because of its entanglement with the function register. To 
better quantify the amount of the quantum correlations created between x and {2/1,2/3}, 



we evaluate the linear entropy el of the qubit xq as [28[: El = 2(1 — Trp^), where the 



factor 2 stems out from the normalization condition and pi is the square of the reduced 
density matrix P{ Xo>x > }(y, y) of the electron of the qubit xo integrated over y. We find that 
El = 0.999 and therefore a maximal correlation between the two registers of the quantum 
circuit is build up, that unambiguously proves the quantum nature of the simulated circuit, 
as required by Shor's algorithm. Once the logical state probabilities of the qubit x are 
known, the latter are combined with the qubit x\ in the zero state and then, as required, the 
order of the argument bits is inverted. This procedure allows one to obtain the binary output 
of the circuit already discussed in Sect. Ill A lj. namely 00 and 10. The first is found with 
a probability of 50,1% and represents the expected failure of the Shor's algorithm, whereas 



the second is obtained with a 49,9% pro 



3ability and leads a to successful determination of 



the order r,. As theoretically expected |3[, failure and success have equal probabilities. 
These outputs indicate an almost ideal performance of the quantum algorithm. 
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FIG. 5: (Color online). Density matrix P{ XOjX > }(y, y) of the electron of qubit xo after the modular 
exponentiation. The diagonal elements describe the density probability of finding the electron in 
the point y along the wire or 1. Here the moduli of P{ XQiX ' j(y, y) are plotted. Note that the 
curves reported in the left panels refer to the left ordinate axes, while the ones reports in the right 
panels refer to the right ordinate axes. 
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B. C=2 



The numerical investigation of the compiled version of Shor's algorithm with C=2 and 
the evaluation of the function ln2[2 a mod 15] in the function register required the evaluation 
of the time evolution of all of the qubits of the registers: Xq, x\, yo, yi- The density 
matrix describing the global system is displayed in Figure [6] at three different time steps. 
The output quantum state describes, with a fidelity of 0.89, the product of two maximally 
entangled Bell pairs, as theoretically expected. The argument register outputs are reported 
in Figure [3, where the reduced density matrix P{ Xo ,x' ,x 1 ,x' 1 }iy 1 y, y', y') of the couple of the 
electrons of qubits {xo,x±} is displayed. The argument register is almost maximally mixed 
as a consequence of the entanglement with the qubits {yo,yi}, as the large value of the 
linear entropy El = 0.976 confirms. The binary output of the algorithm, namely one among 
the possible two-bit responses 00, 01, 10 and 11, is obtained by considering the probabilities 
of the logical state of the qubits {xq,Xi} and then inverting their order. The second and 
the fourth terms yield r=4, which gives correctly the factors 3 and 5 once processed in the 
classical Euclid's algorithm; on the contrary, the first value corresponds to a failure mode 
whereas the third on leads to trivial factors. All of the outcomes have exactly the same 
probability to happen, which consequently means that the routine has a success rate of 50% 
like the previously-discussed case of C—ll. 

Although the binary outputs of the argument registers clearly indicate that the algorithm 
performance is near-ideal both for C—ll and C=2, the efficiency of the quantum networks 




FIG. 6: (Color online). The density matrix of the qubits {xq, x\, yo, yi} evaluated at three different 
time steps: a) input, b) initialization procedure, and c) after the modular exponentiation. Note 
that these matrices have been obtained from the full density matrices of the electrons of the qubits 
{xq, xi, yo, yi} by integrating over the variables defining the position of three carriers along the 
channels. Here the moduli of the density matrix elements are plotted. 
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FIG. 7: (Color online). Density matrix P{ XQ)X ' 0)Xl ^' 1 }{y^ y, y', y') of the couple of electrons of qubit 
{xo,xi} after the modular exponentiation. The diagonal elements describe the density probability 
of finding the two electrons in the points y and y' along the wire or 1. Here the moduli of 

p{ :C o,a ; ' ,»i,4}(y>y>y / )y / ) are plotted. 

is slightly different in the two implementations investigated, as evidenced by the estimated 
fidelity and degree of entanglement between the registers. More specifically, the quantum 
circuit performance gets worse moving from C=ll to C=2. Such a behavior, at a first 
sight, appears to be surprising, due to the larger number of one- and two qubits quantum 
gates numerically simulated in the former implementation. A possible explanation could 
amazingly bring up the small errors inherent the tuning of the quantum gates. In fact, the 
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flaws of the one -and two-logical operations could counterbalance each other resulting as a 
net effect a higher efficiency of quantum circuits. 

V. CONCLUSIONS 

Shor's algorithm highlights the potential power of quantum computation and nowadays 
its realization in scalable structures represents one of the main challenges of quantum in- 
formation science. Only in recent years experimental demonstrations of this algorithm have 
been given in some physical scenarios ranging from NMR to photon qubits. Nevertheless, 
the quantum nature of the processes and/or the scalability of the investigated systems in 
these approaches are questionable. 

In this paper we have introduced and numerically simulated an implementation of the eas- 
iest meaningful example of the Shor's algorithm, that is the factorization of N=15, through 
co- primes C=ll and C=2. The idea we have proposed exploits the coherent SAW-assisted 
transport of electrons in networks of coupled quantum wires, and has a great potential in 
view of its integrability with conventional microelectronics and of its scalability to more 
complex systems containing many qubits. Quantum information is processed by means of a 
sequence of one- and two-qubit gates, materialized by means of an electronic beam splitter 
and phase shifter and a Coulomb coupler, respectively. Their experimental realization in 
semiconductor quantum wires is very challenging since it requires the use of frontier semicon- 
ductor technology. Only in the last years, prototype blocks mimickin g sin gle-qubit rotations 



in a couple of ID channels have been experimentally demonstrated [22|, |29|, |30|. In partic- 
ular, the switching of coupled- quantum wire qubit characteristics has been explored 22]. 
Furthermore, Fischer et al. controlled the coupling between two modes of a couple of ID 
channels, obtained exploiting the two minima of the conduction-band edge in the growth 



direction of a GaAs 2D electron gas 30] or two vertically-coupled 2D electron gases [29]. No 
experimental evidence of two-qubit operation in quantum wires networks has been achieved 
so far. On the other hand, the coherent manipulation of charge states in two spatially 
separated double quantum dots integrated in a GaAs/AlGaAs heterostructure has been re- 



alized 
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321 ] . Specifically, two-qubit operations (swap and controlled-rotation) have been 
successfully implemented. 

We stress that the protocol here proposed for the order-finding routine at the heart of 
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Shor's algorithm represents a"'non standard" implementation of the quantum circuits as 
commonly used in the literature for quantum factorization [l, 4|. Such an implementation 
keeps the basic features of the original algorithm (i.e., "massive parallelism" given by the 
entanglement between the quantum registers and binary output), and also allows for a simple 
network with a lower number of fundamental gates. This makes the numerical simulation of 
the presented protocol less demanding and could also have interesting perspectives on the 
full-scale realization of Shor's algorithm. 

The high efficiency of the quantum processes simulated is shown by the large values 
obtained for fidelities. Furthermore, also the success rate of the algorithm are close to 
its ideal value, in agreement with recent experimental investigations [7]]. The algorithm 
performance is even more noteworthy if we do consider the good but not ideal geometry of 
the logic gates and compare our data with those of the near-ideal case. This behavior is a 
clear signature of the robustness of the algorithm, which is also able to accommodate small, 
but non negligible errors coming from the fabrication and tuning of the quantum gates. The 
capability to take into account small deviations from ideality is certainly a plus, that makes 
the algorithm to compare favourably to any of its experimental implementations. In fact, 
it gives the opportunity to let the device work correctly even in presence of unavoidable 
environmental decoherence effects, always present even at low temperatures. 

Since the recent developments in nanostructure fabrication opened new scenarios in scal- 



able electronic quantum computation 10 
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33|, the promising results here presented in- 



dicate a fruitful guideline for the research in quantum information science. Specifically, this 
work highlights a peculiar physical architecture which could become, in a near future, a 
powerful mean to implement a broader variety of quantum algorithms and therefore to fully 
exploit the whole potential of quantum computation. 
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